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ABSTRACT 

We report on the results of observations in the CO(l-O) transition of a complete 
sample of Southern, intermediate redshift (z = 0.2 - 0.5) Ultra-Luminous Infra-Red 
Galaxies using the Mopra 22m telescope. The eleven ULIRGs with Lpm > 10 12 5 L Q 
south of S — —12° were observed with integration times that varied between 5 and 24 
hours. Four marginal detections were obtained for individual targets in the sample. 
The "stacked" spectrum of the entire sample yields a high significance, 10ct detection 
of the CO(l-O) transition at an average redshift of z = 0.38. The tightest correlation 
of Lfir and Leo f° r published low redshift ULIRG samples (z < 0.2) is obtained 
after normalisation of both these measures to a fixed dust temperature. With this 
normalisation the relationship is linear. The distribution of dust-to-molecular hydrogen 
gas mass displays a systematic increase in dust-to-gas mass with galaxy luminosity 
for low redshift samples but this ratio declines dramatically for intermediate redshift 
ULIRGs down to values comparable to that of the Small Magellanic Cloud. The upper 
envelope to the distribution of ULIRG molecular mass as function of look-back time 
demonstrates a dramatic rise by almost an order of magnitude from the current epoch 
out to 5 Gyr. This increase in maximum ULIRG gas mass with look-back time is even 
more rapid than that of the star formation rate density. 
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1 INTRODUCTION 

The Universe has evolved dramatically with cosmic time. 
The star formation rate density, in particular, first rose in 
the first few Gyr and then declined by more than one order 
of m agnitude in the past 10 Gyr (eg. iHopkins fc Beacoml 
2006). At the same time, the galaxies hosting star formation 
have systematically changed from being dominated by the 
highest mass syst ems in the distant pa st to a more equal mix 
of galaxy masses (|Heavens et al.ll2004l ) in the most recent 0.5 
Gyr, where a representative local volume has been sampled. 
The processes driving this dramatic evolution will only be 
understood if there is a more complete documentation of 
the time evolution of the major baryonic constituents of the 
cosmos. 

An important baryonic constituent to consider in the 
context of galaxies is the molecular gas mass; since this is a 
prerequisite for star formation and as such plays a key role 
in determining its rate and location. While ideally Nature 
would have provided a direct tracer of molecular hydrogen 
mass under all physical conditions, this is unfortunately not 
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the case; making it necessary to infer molecular gas mass 
from indirect measures. An important proxy for molecular 
hydrogen is the detection of carbon monoxide (CO) emis- 
sion lines. Despite the fact that the CO emission lines are 
often highly self-opaque, high resolution studies of individual 
molecular clouds and nearby external galaxies have permit- 
ted at least rough calibration of the total molecular hydro- 
gen mass that is statis tically associated wi th an observed 
CO line luminosity fee. iDickman et ai1 ["l986). The situation 
is complicated by a significant statistical scatter as well as 
the likely dependencies of this association on many other 
factors. A good recent overview of s uch calibration factor s 
for the CO(1-0) transition is given in lTacconi et al.l {2008). 

The relatively high brightness of the CO (1-0) transi- 
tion ha s permitted its rou tine detection in emission at z ~ 
0.1 (eg. lChung erai1l2009l) in a few cases out to redshifts as 
large asz~ 0.4 (jGeach et al.ll2009l .l 201 ll ). and most recently 
in ve ry luminous galaxies at significantly hi gher redshifts as 
well (|Riechers et al.|[2010l : llvison et al.ll201ll ). Higher J tran- 
sitions of CO have been detec ted out to redshifts of 1 < z 
< 3.5 (eg. iGenzel et al.ll2010D . In this paper we address a 
major gap in our observational knowledge of the molecular 
gas mass associated with galaxies by targeting a complete 
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Southern sample of the eleven most luminous FIR detected 
galaxies at redshifts 0.2 < z < 0.5 for deep integrations with 
the Mopra 22m telescope. Only a hand ful of detections have 
yet been made in this r edshift range (|Solomon et al.|[l997l ; 
'beach et al.1 120081 l201ll) . Our study complem ents a simi- 
lar program undertaken bv lCombes et all l|201ll ) directed at 
30 Northern and equatorial targets. Together these studies 
provide some insights into the evolution of the molecular gas 
mass in ULIRGs during this pivotal era in cosmic history. 

The current paper is organised as follows. We begin 
with a brief definition of the sample in describe the ob- 
servations and data reduction methods in Sj3] and present 
the results in A more extensive discussion of the results 
in a cosmological context is deferred to a subsequent pub- 
lication. In this paper we adopt a flat cosmological model 
with Q\ = 0.73 and a Hubble constant of 71 km s~ 1 Mpc _1 
|Hinshaw et alj|2009h . 



2 SAMPLE DEFINITION 

Studies of molecul ar gas in the local Universe (eg. 
ISolomon et "all ll997T l have demonstrated that the largest 
gaseous reservoirs are apparently associated with the galax- 
ies most luminous in the Far Infrared (FIR). This correlation 
provides an opportunity to sample the molecular content at 
higher look-back times by targeting Ultra-Luminous Infra- 
Red Galaxies (ULIRGs, defined by log(LF/i?) > 12) in the 
relevant redshift interval for deep integrations. Estimates of 
the total molecular gas mass can then be made on t he basis 
of the molecular mass function |Keres et al.l 1200*31 ). which 
may evolve, and the space density of the observed target 
population. While this approach is no substitute for com- 
plete sampling of the molecular mass function at each look- 
back time, it provides a first indication of possible trends. 

Our sample was defined by considering all galaxies 
tabulated within the NASA/IPAC Extragalactic Database 
(NEDj] with FIR detections at 60 and 100 /im and spec- 
troscopic redshifts in the range z = 0.2 - 0.5 that exceeded 
a limiting FIR luminosity log(LF/fl) > 12.5. The FIR lu- 
minosity has been defined by Lfir = 4 7rD 2 CC Ffir, 
in terms of the lumin osity distance Dl, a Co lor Correc- 
tion factor, CC = 1.42 ijSanders fc Mirabell[T99l ) and a FIR 
flux density, F FIR = 1.2 6xlQ- 14 (2.58 S 60 +Sioo) W m" 2 
i Sanders fc Mirabel [l99l l. The eleven targets below a Dec- 
lination, S < —12°, that satisfied these requirements are 
listed in Table [1] The Declination cut-off of the sample 
w as chosen to be com plementary to the Northern sample 
of ICombes et~ai] l|20H| j. 



1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



Table 1. The southern ULIRG sample. 



Name 


r> A o-nnn 


DfiConnn 




F00320-3307 


00:34:28.5 


-32:51:13.0 


0.439 


00397-1312 


00:42:15.5 


-12:56:03.0 


0.26172 


00406-3127 


00:43:03.2 


-31:10:49.0 


0.3424 


02262-4110 


02:28:15.2 


-40:57:16.0 


0.49337 


02456-2220 


02:47:51.3 


-22:07:38.0 


0.296 


03538-6432 


03:54:25.2 


-64:23:45.0 


0.3007 


F04565-2615 


04:58:34.7 


-26:11:14.0 


0.490 


07380-2342 


07:40:09.8 


-23:49:58.0 


0.292 


23515-2917 


23:54:06.5 


-29:01:00.0 


0.3349 


F23529-2119 


23:55:33.0 


-21:03:09.0 


0.42856 


F23555-3436 


23:58:06.5 


-34:19:47.0 


0.490 



3 OBSERVATIONS AND REDUCTION 

Observations were carried out on the Mopra 22m telescope"""] 
near Coonabarabran, Australia, between 15 May and 29 Oc- 
tober 2008. Centre frequencies for the target galaxies varied 
between 77 and 91 GHz. The centre frequencies were placed 
in the middle of one of four intermediate frequency bands, 
each of 2 GHz width, and each sampled by 8192 spectral 
channels in two perpendicular polarisations. System temper- 
atures were measured continuously by calibration against a 
noise diode and varied between about 160 and 350 K dur- 
ing useful observing conditions. Calibration of the antenna 
temperature for atmospheric attenuation was updated at 15 
minute intervals by the use of an absorbing paddle at ambi- 
e nt temperature place d over the feed (employing the method 
of lUlich fc Haaslll976h . The telescope pointing was updated 
each hour using the cataloged SiO maser of smallest angular 
separation with the target galaxy, yielding a pointing accu- 
racy of 5 - 10 arcsec, relative to a FWHM beamwidth of 
about 36 arcsec. 

Standard target observing sessions lasted one hour and 
consisted of alternate target and reference spectra, each of 
one minute duration, using a pair of reference locations off- 
set by 15 arcmin to both the East and West of each target 
galaxy. Each target spectrum was calibrated with the aver- 
age of the two adjoining offset reference spectra. Only those 
observing sessions with stable system temperatures below 
about 350 K were retained for further processing. 

Individual one minute integrations on the target were 
subjected to Fourier filtering. Peaks in the Fourier spectrum 
exceeding 10 times its RMS fluctuation level were replaced 
with zero. A more rigorous Fourier clipping was then ap- 
plied to narrow features in the Fourier spectrum. Localised 
Fourier peaks, exceeding the average Fourier amplitude over 
a sliding ten pixel window by 5 times the RMS fluctuation 
level were clipped at this 5a excursion level. Tests of this 
procedure verified that injected signals of the anticipated 
amplitude and linewidth of plausible detections of our tar- 
get galaxies were not significantly degraded by this filtering. 
The typical improvement in the RMS fluctuation level pro- 
vided by this filtering was 50%. The target integrations of 
individual observing sessions were averaged and a second or- 

2 The Mopra radio telescope is part of the Australia Telescope 
which is funded by the Commonwealth of Australia for operation 
as a National Facility managed by CSIRO. 
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Table 2. Source attributes and results. 



Name 


z 




To 


M D 


T 


AT* 


Sco 


Leo 


L'co 


M H2 






(log(L )) 


(K) 


(log(M )) 


(min) 


(mK)t 


(Jy-km/s)t 


(log(L )) 


(log(K- 
km/s pc 2 )) 


(log(M )) 


F00320-3307 


0.439 


12.68 


46 


8.29 


340 


0.63 


<7.9 


<6.58 


<10.89 


<11.10 


00397-1312 


0.26172 


12.67 


52 


8.02 


285 


0.61 


12.6±4.1 


6.32: 


10.63: 


10.79: 


00406-3127 


0.3424 


12.58 


50 


8.03 


283 


0.69 


<7.7 


<6.35 


<10.66 


<10.84 


02262-4110 


0.49337 


12.69 


56 


7.82 


667 


0.60 


<7.8 


<6.69 


<11.00 


<11.13 


02456-2220 


0.296 


12.50 


46 


8.17 


560 


0.51 


<5.8 


<6.10 


<10.41 


<10.63 


03538-6432 


0.3007 


12.58 


49 


8.07 


377 


0.50 


11.1±3.0 


6.39: 


10.70: 


10.89: 


F04565-2615 


0.490 


12.62 


49 


8.11 


915 


0.64 


7.8±3.1 


6.68: 


10.99: 


11.18: 


07380-2342 


0.292 


12.78 


37 


9.05 


1485 


0.40 


<3.8 


<5.90 


<10.21 


<10.52 


23515-2917 


0.3349 


12.54 


46 


8.17 


455 


0.49 


8.5±2.8 


6.37: 


10.68: 


10.89: 


F23529-2119 


0.42856 


12.52 


47 


8.11 


440 


0.48 


<4.1 


<6.28 


<10.59 


<10.80 


F23555-3436 


0.490 


12.67 


16 


8.31 


552 


0.65 


<7.5 


<6.66 


<10.97 


<11.18 


Average 


0.38±0.09 


12.63 


47 


8.29 






4.84±0.48 


6.45 


10.76 


10.96 



^The spectral RMS is listed for a resolution of about 90 km/s. 
* Upper limits are 2a for a 500 km/s linewidth. 



der polynomial baseline was fit to the central 1.75 GHz of the 
band and subtracted from the entire spectrum. All observing 
sessions for individual targets were averaged together with 
an inverse variance weighting. No further baseline fitting of 
any type was applied to the averaged spectrum. 

Total observing times per target varied between 5 and 
24 hours. Spectral smoothing of the final combined spec- 
tra provides a decreasing RMS fluctuation level that scales 
approximately with the square root of bandwidth through 
about 30 km/s. Further spectral smoothing, from 30 to 
90 km/s decreases the fluctuation level by only about 50%, 
rather than the expected 70%, implying that residual sys- 
tematic bandpass effects begin to dominate. The resulting 
RMS sensitivity at a spectral resolution of about 90 km/s 
varied between 0.4 and 0.7 mK in terms of calibrated an- 
tenna temperature, TJ. Main beam brightness temperature 
is given by Tb = TX/rjMB and the main beam efficiency at 
77 - 91 GHz is estimated to be ijmb ~ 0.49 (|Ladd et al.1 
2005). The assumed telescope gain that relates calibrated 
antenna temperature to flux density in this frequency rangeQ 
is 22 Jy/K. 



4 RESULTS 

Final spectra for each of the observed targets are shown 
in Fig. Q] at a spectral resolution of about 90 km/s. Verti- 
cal lines in this figure mark an interval of 1000 km/s cen- 
tred on the reference redshift listed in Table [2] The uncer- 
tainty in the published redshifts is indicated in each panel. 
These can be quite substantial. No high significance detec- 
tions of the redshifted CO(1-0) line were achieved. Upper 
limits to the integrated line-strength can be calculated by 
assuming a representative line-width. Detect ions of CO (1- 
0) in lower redshift samp les of ULIRGs (e.g. ISolomon et al.l 
1 19971 . Ichung et "all I2009T ) suggest that the typical observed 
linewidth is about 300 km/s. We therefore assumed a more 
conservative 500 km/s FWHM for the determination of an 



3 Calibration factors are documented in the Mopra guide 
http://www.narrabri.atnf.csiro.au/mopra/mopragu.pdf 



appropriate noise level. (Flux and mass limits become less 
restrictive the larger the assumed linewidth.) Since we are 
searching for the counterpart of a target with known position 
and velocity, the probability for a false positive detection 
is given by the tail probability of a standard normal dis- 
tribution under the assumption of approximately Gaussian 
noise. A false positive exceeding 2a then has a probability 
of 0.023. Even for our entire sample size of 11 targets, the 
false positive rate at this level of significance is only 0.25, 
which is still significantly less than unity. The 2a upper lim- 
its to integrated linestrength were determined from the mea- 
sured fluctuation level over the entire spectrum (including 
any possible emission signature) after spectral smoothing to 
500 km/s FWHM. The four cases (00397-1312, 03538-6432, 
F04565-2615 and 23515-2917) where the actual line integral 
within the central 1000 km/s of the spectrum exceeds this 2a 
limit are noted as marginal detections in the Table, together 
with their la errors. As noted above, some of the target red- 
shifts have large uncertainties, particularly for F00320-3307 
and F23555-3436. In both these cases, there may be evi- 
dence for CO emission at an offset velocity that may well 
be consistent with the current uncertainties in the systemic 
velocity. 

In the lower right panel of Fig. [1] we present the spec- 
trum obtained by aligning all of the measured spectra to 
relative velocity and forming the simple unweighted aver- 
age. This yields a high signal-to-noise detection of the mean 
CO(1-0) line from our sample, as noted in Tabled 

The corresponding line luminosity is given by, 



47rO|i/ 
c(l + z) 



S(V)dV 



Leo = 1-04 • 10" 



(1+*) 



Sy (L ) 



(1) 



(2) 



for the luminosity distance, Dl , in Mpc, the rest frequency, 
vq, in GHz and the integrated linestrength, 5V, in Jy-km/s. 
Some authors make use of the quantity, 
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Figure 1. Final spectra of CO(l-O) for the target galaxies in the sample shown at a spectral resolution of about 90 km/s. Vertical 
lines mark an interval of 1000 km/s centred on the expected systemic velocity. The published uncertainty in the systemic velocity 
is indicated in each panel. The average stacked spectrum is shown at lower right with the average redshift indicated. 



or 

L'co = 3.26 ■ 10 7 fj- Sv (K - km/s pc 2 ) (4) 
v Q (L + z) 

with Dl, vo and Sv as denned above. For the CO(1-0) tran- 
sition these measures are related simply by, L' co /Lco = 
2.05 • 10 4 (K-km/s pc 2 )/L Q . The utility of the L' co mea- 
sure is that it permits more direct comparison of different 
line transitions since it is formulated as a product of bright - 
ness temperature with surface area (|Solomon et all 1 199 7ft . 
Line emission of the same brightness temperature originat- 
ing from the same region will yield the same L' co , indepen- 
dent of transition or observing frequency. 

Calculation of an associated total hydrogen mass re- 



quires an assumed relationship between this quantity and 
the CO (1-0) line luminosity. A goo d compilation of suc h 
conversion factors is given in Fig. 10 of lTacconi et alj {2008). 
There is clearly a large scatter in the derived conversion 
factors and there may well be underlying dependancies on 
many physical factors, including metallicity, gas surface den- 
sity and particularly gas excitation temperature. A plausi- 
ble conversion factor appropriate for ULIRGs is likely to be 
about 1 Mq per (K-km/s pc 2 ) or about 2 • 10 4 Mq per Lq. 

We plot the distribution of CO(1-0) and FIR luminos- 
ity in Fig. [21 i nclud ing both our own results and those of 
ICombes et al.l (|201 lh as well as previously p ublished studies 
of less distant and less luminous systems bv lSolomon et all 
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Figure 2. The relationship between CO(l-O) and FIR luminosity. 
Our results are indicated with the filled circles and 2a upper limit 
symbols. The large filled pentagon represents our "stacked" result. 
A power law of slope 1.37±0.06 represents a least squares fit to 
the low redshift samples and has a correlation coefficient 0.891. 



jl997h . ICh.ung et all l|2009l ) and lGao fc SolomorJ (|2004h . all 
of which have been undertaken in the CO (1-0) transition. 
Although the low redshift data lie on a well-defined locus in 
the plot, the interme diate redshift data fr om both the cur- 
rent study and that of lCombes et al.l ()201 lh show divergence 
from this simple correlation. A subset of the detections has a 
significant excess CO luminosity, while another subset con- 
sisting primarily of upper limits in CO, is significantly un- 
derluminous. The difference in CO luminosity between these 
two populations is about one order of magnitude. A power- 
law with slope of 1.37±0.06 has been fit by least squares 
to the three low redshift samples and is overlaid in the fig- 
ure. The correlation coefficient of the data points is 0.891. 
As noted by previous authors (eg. IChung et all [2009). this 
relationship is steeper than linear. A single power-law slope 
provides a reasonable representation of this correlation for 
FIR luminosities 9.3 < \og(LFm) < 12.2. As noted pre- 
viously, the associated gas mass for a given CO (1-0) line 
luminosity is likely to vary with the gas kinetic and excita- 
tion temperatures, while Lfir is expected to be very sensi- 
tive to the dust te mperature, Tb, possibly varying as T^~ 6 
l|Soifer et al.lll989l). 



Gao fc Solomon! (|2004T ) have demonstrated a higher de- 



gree of correlation of Lfir with L' co once a correction 
of the FIR luminosity for a varying dust temperature has 
been accounted for. The dust temperature can be estimated 
from the FIR photometry if an emmissivity law is assumed. 
This is often taken to be that of a Planck function multi- 
plied by , fo r frequ ency, v, and power law exponent, /3. 
iLisenfeld et al.l (|2000h have found that the spectral energy 
distributions (SEDs) in their sample of FIR luminous galax- 
ies could be fit with 0.85 < /3 < 1.9. The dust tempera- 
ture, in the Wien approximation, can be estimated from the 



ratio, Rs = S1/S2, of flux densities at two fixed observing 
frequencies, v\ and v-z from, 



(5) 



where R v — v\/u2- For the 60 and 100 fim IRAS bands 
this yields, 

41.7(1 + 0) 



(K) 



(6) 



0.22(/3 + 3) - \og(Rs) 

Since the integrated FIR luminosity is expected to scale as 
T^ +4 , we have considered normalisations of the FIR lumi- 
nosity down to a reference temperature of 25 K using the 
form, 



-((3+4) 



(7) 



The case has also been made that the CO(1-0) line 
luminosity would va ry approximately linearly with Tb 
l|Solomon et al.lll997h if the dust and gas temperatures are 
reasonably coupled, suggesting a similar normalisation of 
the form, 



L'a 



To 
25 



(8) 



We find empirically,that the values (a,/3) — (1,1) minimise 
the dispersion in the temperature- normalised CO (1-0) ver- 
sus FIR luminosity relationship. The resulting relation is il- 
lustrated in Fig. [3] where a power-law of slope 0.99±0.04 has 
been fit by least squares to the three low redshift samples. 
The correlation coefficient with this choice of (a, /3) is 0.902. 
A very similar relation, with a slope 1.03±0.05, is obtained 
for (a, P) = (1,1.5), which may be more representative of the 
FIR SEDs jLisenfeld et al.ll2000l) . The distinction between 
the CO-luminous and CO-poor populations is undiminished 
in the temperature corrected relationship. 

We can estimate the dust mass associated with our tar- 
gets from, 



SuoD L 



{l + z)K vz B uz (Tn) (9) 

where SVo is the observed flux density in an FIR band, Dl 
is the luminosity distance, k uz is the absorption cross sec- 
tion per unit dust mass at the target rest-frame frequency 
v z = ^o(l + z) and B uz (Td) is the Plan ck function at th is 
frequency for a dust temperature, To l)Hildebrandl 1 19831 ) . 
The dust absorption cross section can be approximated by, 



ka = 27.2 



Too" 



(l + 0.625[log(A Mm ) - 2.15] 2 ) (10) 



in units of cm 2 /g for wavelengths expressed in /im between 
40-100 (im. This analytic form fits the tabulated data of 
iDraind (120031 ) to better than 5% over the indicated wave- 
length range. We calculate the dust masses for our sample 
and the comparison samples noted previously using the ob- 
served IRAS 100 /im fluxes and the dust temperatures cal- 
culated with eqn. [6] and /3 = 1.5, listing these in Table [2] 

Once the CO luminosity has been temperature- 
corrected, it seems more likely that a single conversion factor 
to total molecular hydrogen mass, appropria te to the 25 K 
reference temperature, might apply. From iTacconi et al.l 
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Figure 3. The relationship between temperature normalised 
CO(l— 0) and FIR luminosity. Lqq ls normalised by (Tjj/25) 1 and 
Lpm by (Trj/25) , for a reference dust temperature, To = 25 K. 
A power law of slope 0.99±0.04 represents a least squares fit to 
the low redshift samples and has a correlation coefficient 0.902. A 
similar fit is achieved for Leo normalised by (Tjj/25) 1 and Lfir 
by (T D /25) 5 ' 5 . 
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Figure 4. The relationship between molecular hydrogen mass 
and dust mass. A power law of slope 1.09±0.05 is shown as a 
dashed line and represents a least squares fit to the low redshift 
samples and has a correlation coefficient 0.893. Dotted lines la- 
belled with sub-SMC, LMC and MW represent increasing dust- 
to-molecular hydrogren mass ratios of 0.001, 0.003 and 0.01 re- 
spectively. 



(2008) this might correspond to about 3 Mq per (K-km/s 
pc 2 ) or about 6- 10 4 Mq per Lq. Even so, there may well be 
other systematic dependancies and a large intrinsic scatter, 
that add significant uncertainty to this conversion factor. 
The relationship between molecular hydrogen mass and dust 
mass is shown in Fig. [3] A power-law fit to the low redshift 
samples has slope 1.09±0.05 and a correlation coefficient of 
0.893 and is plotted in the Figure. Also shown are diagonal 
lines corresponding to fixed dust to hydrogen mass ratios 
of 0.001, 0.003 and 0.01, which are labelled sub-SMC, LMC 
and MW respectively (cf. lDraine et al.ll2007T ) . The sub-SMC 
designation is used for the lowest curve since the Small Mag- 
ellanic Cloud is estimated to have a mass ratio of 0.002 (eg. 
IWeingartner fc Draine|[200ll ). The correlated locus in Fig. [4] 
defined by the low redshift samples systematically increases 
in dust-to-gas mass ratio with increasing gas mass, from 
SMC-like values at the low end, through Large Magellanic 
Cloud and up to Milky Way values for Mh 2 > 10 10 M Q . 

What is striking in Fig. [4] is that all of th e intermediate 
redshift detections in both the current and IComb es et al.l 
(|20 1 if ) studies depart from the low redshift trend in the 
direction of extremely low dust-to-molecular gas mass ra- 
tio. Our well-defined sample average point of 0.0021±0.0002 
is comparable to the dust-to-hydrogen mass ratio of the 
SMC. In contrast, there are several upper limits in the fig- 
ure which appear to be dust-rich (or gas-poor) by about 
an order of magnitude. The most extreme point from the 
current sampl e is th e source IRAS 0738-2342, while in the 
ICombes et ail (|201lh sample it is IRAS 19104+8436. Two 
addtional sources in the Northern sample, [HB89] 1821+643 
and F 00415-0737, are also discrepant but less extreme. 
While very little is currently known about F 00415-0737, 
iRuiz et all (|2010T ) have recently presented SED fits for IRAS 
0738-2342 and [HB89] 1821+643. They estimate relative 
AGN:Starburst contributions to the bolometric luminosity 
of about 50:50 and 80:20 for these two sources. Both IRAS 
19104+8436 and [HB89] 1821+643 are classified as Seyfert 
1/QSOs, suggesting that L FIR for IRAS 19104+8436 may 
also be dominated by non-thermal rather than dust emis- 
sion. 

It seems plausible that the apparent "gas poor" sub- 
population seen at intermediate redshift in Figs. [2] [3] and [4] 
may simply be a manifestation of AGN contamination, with 
the star formation dominated population demonstrating a 
well-defined trend relative to the low redshift samples, con- 
sistent with a dramatic decline in the dust-to-molecular gass 
mass ratio. 

We consider the distribution of ULIRG molecular hy- 
drogen mass with look-back time, r, in Fig. [5] Since ULIRGs 
are identified by their extreme luminosity in the FIR band, 
where self-opacity and foreground extinction effects are min- 
imal, they represent a population that is easily recognised 
in unbiased surveys of the sky. Both our own study, as well 
as the others shown in the Figure, were drawn from the 
most luminous soures in the IRAS all-sky survey. As such 
they are not expected to suffer from incompleteness at high 
luminosity and despite the fact that they are a rare pop- 
ulation, they permit a useful assessment of the associated 
molecular gas mass that is required to feed the ULIRG phe- 
nomenon. Although the lowest detected gas masses at each 
look-back time simply reflect the sensitivity of the observa- 
tion and the cutoff Lfir of the sample, the upper envelope 
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Figure 5. The relationship between look-back time and ULIRG 
molecular hydrogen mass. An indicative upper envelope to molec- 
ular gas mass, log(Mj^ 2 ) = 10.2 + 0.18r, is plotted demonstrating 
the order of magnitude decline in maximum ULIRG gas mass over 
the past 5 Gyr. 



to mass in Fig. [5] is a significant attribute of the population, 
at least once the sample extends over a representative vol- 
ume of about 10 7 Mpc 3 ; since the most luminous sources 
are clearly the easiest to detect at any redshift. This condi- 
tion should be met for redshifts greater than about 0.03, or 
look-back times exceeding 0.4 Gyr. This local volume satu- 
ration effect is apparent in the Figure for r < 0.5 Gyr. The 
indicative curve drawn in the Figure, 



log(Afo 2 ) = 10.2 + 0.18r 



(11) 



demonstrates the order of magnitude decline in maximum 
molecular gas mass of the ULIRG population during the 
past 5 Gyr. Such a decline is comparable to, but even more 
dramatic, than that seen in the star formation rate density 
iris fc Beacomll2006h which can be well described by, 

log (p.) = -1.76 ± 0.02 + 0.131 ± 0.004r (12) 

for a star formation rate density p* in M©yr -1 Mpc~ 3 for 
r < 8 Gyr. This contrast suggests significant changes in the 
molecular gas mass function at even these modest look-back 
times. 

We will present further analysis of the time evolution 
of molecular gas and other baryonic constituents of galaxies 
in a subsequent publication. 
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